Chirally symmetric but confining dense and cold matter. 
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The folklore tradition about the QCD phase diagram is that at the chiral restoration phase 
transition at finite density hadrons are deconfined and there appears the quark matter. We address 
this question within the only known exactly solvable confining and chirally symmetric model. It 
is postulated within this model that there exists linear Coulomb-like confining interaction. The 
chiral symmetry breaking and the quark Green function are obtained from the Schwinger-Dyson 
(gap) equation while the color-singlet meson spectrum results from the Bethe-Salpeter equation. 
We solve this model at T=0 and finite chemical potential and obtain a clear chiral restoration 
phase transition at the critical value ficr- Below this value the spectrum is similar to the previously 
obtained one at jj. — 0. At fi > the quarks are still confined and the physical spectrum consists 
of bound states which are arranged into a complete set of exact chiral multiplets. This explicitly 
demonstrates that a chirally symmetric matter consisting of confined but chirally symmetric hadrons 
at finite chemical potential is also possible in QCD. If so, there must be nontrivial implications for 
astrophysics. 

PACS numbers: 11.30.Rd, 25.75.Nq, 12. 38. Aw 
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INTRODUCTION 

The structure of the QCD phase diagram is one of the 
keys for our understanding of QCD, strongly interact- 
ing matter and has direct astrophysical implications. It 
is a subject of a significant effort within the heavy ion 
collision programs at different labs. We know that in 
the vacuum as well as at low temperature and density 
QCD is realised in the confining hadronic phase with 
broken chiral symmetry. At high temperature and low 
density the system is deconfined and chiral symmetry 
is restored. Above the critical temperature one experi- 
mentally observes a strongly interacting plasma with dy- 
namics that is not yet well understood. Only at the 
very high temperatures one possibly approaches a per- 
turbative regime of the weakly interacting quark-gluon 
plasma. Just in the opposite limit of the asymptotically 
large densities and small temperature one believes that 
the perturbative giuonic interaction is also adequate and 
consequently one has, depending on the number of active 
flavors the color superconducting phase (the color-flavor 
locking phase with TV/ = Nc). 

There is a folklore tradition that upon increasing tem- 
perature and density from zero the deconfining and chi- 
ral restoration transitions coincide and beyond the semi- 
circle in the T — fj, plane one obtains a deconfining and 
chirally restored matter. From the theoretical side we 
know from the 't Hooft anomaly matching conditions 
[H that at zero temperature and density in the confin- 
ing phase chiral symmetry must be necessarily broken 
in the vacuum. The converse is not required, of course. 
Whether the chiral restoration and deconfinement tran- 
sitions coincide at finite temperature and zero chemical 
potential can be answered on the lattice. It was believed 
that the same temperature triggers both transitions 



though recently there have appeared indications that it 
could be not so ^3]. Unfortunately lattice cannot help 
us at finite chemical potential due to the notorious sign 
problem. Hence one has to rely on models of QCD. Typ- 
ically models that give us information about the phase 
structure are of the Nambu and Jona-Lasinio type, for 
an overview and references see [4|. An obvious draw- 
back of all these models is that they are not confining. 
Then there is no basis to conclude that above the chiral 
restoration point at finite chemical potential one obtains 
a chirally symmetric and deconfining quark matter. 

Recently McLerran and Pisarski presented qualitative 
large Nc arguments showing that at a reasonably large 
chemical potential and low temperature there might exist 
a confining but chirally symmetric phase This sug- 
gestion would dramatically change the existing paradigm 
and requires urgent efforts to verify whether it could be 
indeed the case. This suggestion is in conflict with the 
naive intuition that once the hadrons are confined chi- 
ral symmetry should be necessarily broken. There are 
no such limitations from the 't Hooft anomaly matching 
conditions at finite density, however, and such an intu- 
ition is in fact an artefact of existing simple models of 
chiral symmetry breaking and confinement 0, 01 • 

Clearly we cannot solve QCD and answer this ques- 
tion from first principles. What can be done is to verify 
this conjecture within a solvable manifestly confining and 
chirally symmetric model. Such a model could provide 
insight into physics that would drive such an unexpected 
situation. If confirmed within such a model, this sce- 
nario could also be realised in QCD and further theo- 
retical and experimental efforts to clarify this important 
question would be required. 

There exists only one known manifestly chirally- 
symmetric and confining model in four dimensions that 
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is solvable sometimes called Generalized Nambu and 
Jona-Lasinio model (GNJL). This model can be consid- 
ered as a generalization of the 1+1 dimensional 't Hooft 
model, that is QCD in the large limit 0. Once the 
gauge is properly chosen in 1+1 dimensions the Coulomb 
interaction becomes a linear confining potential. Then 
this potential properly represents gluonic degrees of free- 
dom in 1+1 dimensions. In four dimensions there are dif- 
ferent gluonic interactions and it is hopeless to solve even 
large Nc QCD with full gluodynamics. It is postulated 
within the GNJL model that there exists a linear confin- 
ing potential of the Coulomb type in four dimensions [1^ . 
This model represents a simplification of large Nc QCD. 
The model is exactly solvable. The chiral symmetry 
breaking and the properties of the Goldstone bosons can 
be obtained from the solution of the Schwinger-Dyson 
and Bethe-Salpeter equations [13, El [H, [11 Q [H. 



The cornplete spectrum of qq mesons has been obtained 
in refs. 3, 



17| , which exhibits restoration of chiral sym- 
metry in excited hadrons, for a review see [l8| . 

It is instructive to review the chiral symmetry and con- 
fining properties of quarks in this model, which are sim- 
ilar to those in the 't Hooft model. Upon solving the 
Schwinger-Dyson equation for the quark Green functions 
one obtains chiral symmetry breaking in the vacuum, 
i.e. a nonzero quark condensate and a dynamical mass 
M{q) for the quarks, that is strongly momentum depen- 
dent and vanishes at large momenta. This single quark 
Green function is infrared-divergent, which means that 
the single quark cannot be observed and hence is con- 
fined. The infrared divergence is canceled exactly in the 
Bethe-Salpeter equation for the observable color-singlet 
hadrons and the physical spectrum consists only of color- 
singlet hadrons. This model allows for a solution at finite 
chemical potential and therefore can be used as a tool 
to check whether it is possible or not to have confining 
but chirally symmetric matter at finite chemical poten- 
tial. Below we solve it and obtain a complete spectrum 
of mesons at zero temperature and finite quark chemical 
potential. We demonstrate that below the critical chem- 
ical potential, the spectrum is similar to the one 
previously obtained at zero chemical potential [3, [l3] ■ 
At fi = ficr the chiral restoration phase transition hap- 
pens. At /i > /icr hadrons are still confined and the 
spectrum represents a complete set of chiral multiplets. 
This explicitly demonstrates that a chirally symmetric 
but confining phase exists within this model and opens a 
perspective to study such a possibility in QCD. 



CHIRAL SYMMETRY BREAKING AND 
BETHE-SALPETER EQUATION FOR MESONS 
IN VACUUM 

Throughout this paper we work in the chiral limit. The 
two flavor version of the model is considered. The global 



chiral symmetry of the model is U{2)l x U{2)ii, because 
in the large Nc world the axial anomaly is absent. The 
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model is described in great detail in references [17 
here we give only an overview of this model in vacuum. 
It is postulated within the model that there is only a 
linear instantaneous Lorentz-vector interquark potential 
that has a Lorentz structure of the Coulomb potential, 
i.e., it is a density-density interaction: 



g^.o9.oS'''V{\x 



2/1); -^V{r)=crr, 
(1) 

where a, b are color indices. The Fourier transform of this 
potential and any loop integral are infrared-divergent. 
Hence infrared regularization is required and any physi- 
cal observable, such as a hadron mass, must be indepen- 
dent of the infrared regulator ///^ in the infrared limit 
(i.e. when this regulator approaches 0, /i/fl — > ). There 
are several physically equivalent ways to perform this in- 
frared regularization. Here we follow Ref. |llj] and define 
the potential in momentum space as 



V{p) = 



8na 



2 ■ 
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Then this potential in the configuration space con- 
tains the required or term, the infrared-divergent term 
—(7/(1111 well as terms that vanish in the infrared limit. 
The Dirac operator for the dressed quark is 

D{pq,p) = iS^'^{po,p) = Do{po,p) - T:{po,p), (3) 

where Dq is the bare Dirac operator. Parametrising the 
self-energy operator in the form 



^ip)^Ap + ii-p)[Bp~p], 



(4) 



where functions Ap and Bp are yet to be found, the 
Schwinger-Dyson equation for the self-energy operator 
in the rainbow approximation,which is valid in the large 
Nc limit for the instantaneous interaction, is reduced to 
the nonlinear gap equation for the chiral angle ipp, 



Ap cos ipp — Bp sin ipp — 0, 



where 



1 f d^k 



(5) 



(6) 



1 f (Pk ' ^ 
Bp ^ P+2J iP-k)V (p-k) cos ^k- (7) 

The functions Ap,Bp, i.e. the quark self-energy, are 
divergent in the infrared limit. 
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Bp — 



2AtiR 



cos ipp 



+ 

' ^p J 



(8) 



(9) 



where A^^ and are infrared- finite functions. This im- 
plies that the single quark cannot be observed and the 
system is confined. However, the infrared divergence can- 
cels exactly in the gap equation ^ so this equation can 
be solved directly in the infrared limit [26*]. The chiral 
symmetry breaking is signalled by the nonzero quark 
condensate and by the dynamical momentum-dependent 
mass of quarks 



Nc 

{qq) = -— 



dp 



simpp, M{p) — ptaiKpp. 

(10) 

The dynamical mass is finite at small momenta and van- 
ishes at large momenta. Both these quantities were first 
obtained in ref. and repeat edly reconfirmed in all 
subsequent works on this model [ij, [H, [H, Q M, M, 
17[. The numerical value of the quark condensate is 
(qq) = (-0.2317^)3. 

The homogeneous Bethe-Salpeter equation for a quark- 
antiquark bound state in the rest frame with the instan- 
taneous interaction is 



d-^q 



y{\p- q\) 7oS{qo + 711/2, p- q) 

(11) 



(2^)4 

X x{m,q)S{qo - ■m/2,p- q)'^Q 



Here m is the meson mass and p is the relative mo- 
mentum. The Bethe-Salpeter equation can be solved by 
means of expansion of the vertex function x{'>^tP) into 
a set of all possible independent amplitudes consistent 
with I,J^^. Then the Bethe-Salpeter equation trans- 
forms into a system of coupled equations 17|. The in- 



frared divergence cancels exactly in these equations and 
they can be solved either in the infrared limit or for very 
small values of the infrared regulator. 

In the limit of vanishing dynamical quark mass M{p) — 
(i.e., ipp — 0), the Bethe-Salpeter equation transforms 
into systems of coupled equations that exactly fall into a 
complete set of chiral multiplets p^ : 

J = 



(0,0) 
(l/2,l/2)a 
(1/2,1/2)6 
(0,1)0(1,0) 



0, J— 

1, J-+ 
1,J++ 
1,J++ 



0, J++ 

0, J++ 
G,J-+ 

1, J"" 



(13) 



J = 2k-l, k=l,2,. 



(0,0) 


0, J++ ^ 




(l/2,l/2)a 


1,J+-- 




(1/2,1/2)6 


1,J— ^ 


^ 0, j+- 


(0,1)0(1,0) 







(14) 

Here /, J^^ is a standard notation for isospin /, spin 
J as well as for both spatial and charge parities P and 
C. The sign < — > indicates that both given states belong 
to the same representation and that the Bethe-Salpeter 
equations for these states are identical. 

The complete spectrum of th e qq mesons in the vac- 
uum can be found in refs. [11, 17 1 and exhibits a fast 
chiral restoration with increasing J. A reason for this 
restoration is that with larger J the radial wave func- 
tion vanishes at small relative momenta, and hence the 
large chiral symmetry breaking dynamical mass M{p) 
becomes irrelevant and asymptotically at J — * oo the 
Bethe-Salpeter equation exactly decouples into a set of 
equations (fT3l) - (fTil) . 



INCLUSION OF A FINITE QUARK CHEMICAL 
POTENTIAL AT ZERO TEMPERATURE 

The gap equation 

Now we are in a position to include into this model 
a finite quark chemical potential at zero temperature. 
We can straightforwardly do it only at zero temperature 
because to employ thermodynamics at finite temperature 
one would need to know a temperature dependence of the 
confining interaction. In addition there is no screening 
of the confining interaction at zero temperature in the 
large Nc limit since in this limit there are no vacuum 
quark loops and no the particle-hole modifications of the 
gluonic propagator due to Debye screening Q. 

The effect of a finite quark chemical potential at zero 
temperature reduces to the Pauli blocking of all occupied 
levels. Denoting the Fermi momentum of quarks as pf 
one has to replace the vacuum density matrix v{p)v''{p) 
by the density matrix in the medium with all quark levels 
occupied up to the Fermi momentum: 



(l/2,l/2)a : l,0-+^0,0^ 
(l/2,l/2)b : 1,0++ ^0,0- 

J = 2k, k=l,2,... 



(12) 



p{p) = Q{pf - p)u{p)v){p) +v{p)v''{p). (15) 

Hence one has to remove from the integration both in the 
Schwinger-Dyson (gap) and Bethe-Salpeter equations all 
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FIG. 1: Quark condensate in units of a'^^'^ as a function of 
the Fermi momentum, which is units of \/a. 

quark momenta below since they are Pauli-blocked. 
The modified gap equation at any p > is then the same 
as in ([5]) - (O, but the integration starts not from fc = 0, 
but from k = pf. The points p < Pf are irrelevant for the 
hadron which is always on top of the Fermi sea. Similar, 
in the expression (jTU]) for the quark condensate the lower 
limit of the integration is nowp/, rather than 0. As in the 
vacuum, we are in the position to solve numerically the 
nonlinear gap equation at arbitrary pf, study the chiral 
restoration phase transition at the critical pY , and, given 
a dressed quark Green function, solve the Bethe-Salpeter 
equation for the spectrum below and above the critical 
Fermi momentum. 

In Fig. 1 we show the quark condensate (qq) as a 
function of p/. One observes an obvious phase transi- 
tion at the critical Fermi momentum pY — 0.109^/ct. 
This is similar to what has been observed within the 't 
Hooft model The reason for this chiral symmetry 

restoration phase transition is the same as in the stan- 
dard Nambu and Jona-Lasinio model [i^] or in the 't 
Hooft model [l9| : Once a considerable part of the quark 
self-energy interaction is removed by the Pauli blocking, 
then there is not enough strength in the gap equation 
to generate a nontrivial solution with the broken chiral 
symmetry in the vacuum. 

In order to see this physics explicitly consider solutions 
of the gap equation for the chiral angle tpp at different 
fixed Fermi momenta p/, that are shown in Fig. 2. One 
obviously observes a universal approaching of the non- 
trivial solution of the gap equation ipp ^ to its trivial 
value ifp — once the Fermi momentum pf approaches 
a critical value pY- Note that the infrared point in the 

gap equation ([5])- ([7]) is k — p (i.e., where the argument of 
the potential approaches 0). It directly follows from the 
self-energy integrals (I6|)-(l7|). For any p > pf this point 
is always within the integration interval from k — pf to 
k — oo. However, while the linear confining potential is 
infrared-singular itself, the gap equation ^ is not, be- 



FIG. 2: Chiral angle (pp as a function of the momentum p at 
different fixed values of the Fermi momentum pf below the 
critical value. Both momenta are units of ^/a. 

cause the infrared divergence cancels out exactly, as it 
follows from ©-((HI). Hence the contribution to the gap 
equation comes not only from the vicinity of the infrared 
point p = k, but actually from the whole integration in- 
terval. Once a critical part of this integration interval is 
removed by the Pauli blocking, the nontrivial chiral sym- 
metry breaking solution of the gap equation abruptly dis- 
appears. A Pf dependence of the dynamical mass M(ji) 
is shown in Fig. 3. 

Summarizing, we stress that the symmetry gets re- 
stored in our case not because the potential gets screened. 
In the large Nc limit it cannot be screened and the 
string tension a remains constant at any pf. It gets re- 
stored exclusively due to the Pauli blocking that reduces 
a "strength" of the gap equation below the critical value. 
At each fixed value of a there always exists such a criti- 
cal reduction of the integration interval due to the Pauli 
blocking so that the nontrivial chiral symmetry breaking 
solution vanishes. This directly follows from the fact that 
all dimensional quantities in our task are given in units 
of cr and all results are valid for any a, large or small. 

Above the critical value pY there is no nontrivial chi- 
ral symmetry breaking solution and the only solution is 
trivial, (pp = 0. This follows from the numerical inte- 
gration of the highly nonlinear gap equation. It would 
be difficult, if impossible, to obtain this critical Fermi 
momentum analytically. 

Hence the chiral symmetry gets restored, ipp = 0, with 
vanishing quark condensate, (qq) — 0, and dynamical 
mass of quarks, M{q) = 0, as it follows from pil)) . At 
(pk — the Lorentz-scalar self-energy of quarks vanishes 
identically, Ap = 0, see eq. (6). 

What crucially distinguishes this model from the NJL 
model is that the quark is still confined, even in the chi- 
rally restored phase. This is because for any p > pf the 
Lorentz spatial-vector self-energy integral Bp docs not 
vanish at ipk = (see eq. ([7])) and is in fact infrared- 
divergent, see Hence the single quark is removed 
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FIG. 3: Dynamical mass of quarks M(p) as a function of the 
momentum p at different fixed values of the Fermi momentum 
Pf below the critical value. Both momenta are units of ^/a. 

from the spectrum at any chemical potential. 



Meson spectrum 

This infrared divergence of the single quark Green 
function cancels exactly, however, in the color-singlet 
quark-antiquark system [iTj and the bound state mesons 
are finite and well defined quantities. Like in the gap 
equation, the only modification of the Bethe-Salpeter 
equation is that the integration in q starts not from 0, 
but from q — pf. 

The complete spectrum of qq mesons at the Fermi mo- 
mentum Pf essentially below the critical value, pf — 
O-OS-^/ct, is shown on Fig. 4. This spectrum is similar 
to the previously obtained one in the vacuum in refs. 
[liL [TtI I . Four Goldstone bosons with the quantum num- 
bers I = 1,0 ^ and / = 0,0 ^ are well seen. We re- 
mind that there are no vacuum fermion loops in this 
large Nc model and hence the U{1)a symmetry is broken 
only spontaneously. The spectrum exhibits approximate 
restoration of the chiral symmetry in excited hadrons, for 
details we refer to [l3l and for a review to ref. [l8| . 

The spectrum above the critical Fermi momentum, i.e. 
for Pf — is shown in Fig. 5. This spectrum 

is qualitatively different as compared to Fig. 4. All the 
states are in exact chiral multiplets. A fundamental sym- 
metry reason of this degeneracy is rather obvious: When 
the vacuum is chiral- invariant, i.e., we are in the Wigner- 
Weyl mode, the bound state spectrum (if it exists) can 
only be in exact chiral multiplets (2l| . In our case in 
the Wigner-Weyl mode the quarks are confined and the 
bound states do exist. Consequently the resulting color- 
singlet spectrum is manifestly chirally-symmetric. This 
also directly follows from the analytical symmetry prop- 
erties of the Bethe-Salpeter equation in the chiral sym- 
metry limit M(p) = (l7l|. 

Though the chiral symmetry is manifestly restored. 



one observes finite- energy well defined hadrons. Obvi- 
ously the mass generation mechanism in these hadrons 
has nothing to do with the chiral symmetry breaking and 
is not related with the quark condensate. The mass gen- 
eration mechanism for these chirally symmetric hadrons 
comes from the manifestly chirally symmetric dynamics 
and is similar to the high-lying states in the chiral sym- 
metry broken phase on Fig. 4. 

On Fig. 6 we show an evolution of the ground state 
pion and the sigma mesons as a function of the Fermi 
momentum. 

CONCLUSIONS 

As a conclusion we have demonstrated that it is pos- 
sible to have at finite chemical potential a confining but 
chirally symmetric matter consisting of chirally symmet- 
ric hadrons. For that we used the only known exactly 
solvable model that is confining and manifestly chirally 
symmetric. Certainly this model is not QCD and perhaps 
confinement mechanism in QCD is more complicated. 
However, what we have shown, is "a matter of princi- 
ple demonstration". We have also clarified perhaps a 
generic mechanism of the phenomenon. Namely, the chi- 
ral symmetry restoration means that the Lorentz-scalar 
part of the quark self-energy vanishes. But it does not 
require yet that the other parts of the quark self-energy 
vanish either. If they do not and the quark is still off- 
shell, then it cannot be observed and such a single quark 
is confined. At the same time the color-singlet hadrons 
are well-defined and finite on-shell systems. 

Whether this happens in QCD or not is still an open 
question but it definitely suggests that there are no rea- 
sons to believe that deconfinement and chiral restoration 
necessarily coincide at finite density. If such a chirally 
symmetric but confining phase does exist in QCD, then 
it will imply dramatic modifications of the QCD phase 
diagram. Schematically this diagram would look like on 
Fig. 7. 

It will also have significant implications for astro- 
physics: The interactions of these chirally-symmetric 
hadrons can be only of short-range as they decouple from 
the Goldstone bosons and their weak decayrate is quite 
different since their axial charge vanishes |22| . 
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